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We report the results of a weakly interacting massive particle (WIMP) dark matter search using 
the full 80.1 live-day exposure of the first stage of the PandaX experiment (PandaX-I) located in the 
China Jin-Ping Underground Laboratory. The PandaX-I detector has been optimized for detecting 
low-mass WIMPs, achieving a photon detection efficiency of 9.6%. With a fiducial liquid xenon 
target mass of 54.0 kg, no significant excess events were found above the expected background. 

A profile likelihood ratio analysis confirms our earlier finding that the PandaX-I data disfavor all 
positive low-mass WIMP signals reported in the literature under standard assumptions. A stringent 
bound on a low mass WIMP is set at WIMP mass below 10GeV/c 2 , demonstrating that liquid 
xenon detectors can be competitive for low-mass WIMP searches. 

PACS numbers: 95.35.+d, 29.40.-n, 95.55.Vj 


INTRODUCTION 

The existence of gravitationally attractive “dark mat¬ 
ter” that dominates the matter composition of the uni¬ 
verse has been firmly established based on overwhelming 
evidence from astronomical and cosmological observa¬ 
tions if]. Whether such abundant matter consists of yet 
unknown elementary particles remains one of the most 
pressing scientific questions. There are strong theoret¬ 
ical motivations for the existence of beyond the Stan¬ 
dard Model physics, many of which naturally predict 
new stable neutral particles at the electroweak symmetry 
breaking scale with weak interactions, generically named 
weakly interacting massive particles (WIMPs) [5] [3]. 
WIMPs are a leading dark matter (DM) candidate where 
weak interactions between WIMPs and ordinary mat¬ 
ter allow for a direct search for these particles through 
particle physics experiments. In recent decades, direct 
searches of WIMP interactions with terrestrial detectors 
have been carried out in deep underground laboratories 


worldwide with ever increasing discovery power [4]. 

Since 2008, a number of underground direct detec¬ 
tion experiments have reported signals that could be 
interpreted as WIMP interactions within the detector. 
Among those are the DAMA/LIBRA experiment using 
Nal(Tl) crystals [5], the CoGeNT experiment j5] us¬ 
ing point-contact Ge detectors, the CRESST-II experi¬ 
ment [7\ using cryogenic CaWO .4 bolometers (excess not 
reproduced in the recent experiment 0 ). as well as the 
CDMS-Si experiment using cryogenic Si bolometers [5]. 
Although the claimed signals are not generically con¬ 
sistent, they all point to low to median WIMP mass 
in the range of 10 to 50GeV/c 2 . On the other hand, 
the ZEPLIN-III [TO], XENON-100 [II], LUX [TO], and 
PandaX-I m experiments utilizing xenon, the DarkSide- 
50 experiment using argon [TO] , the SuperCDMS mm 
and CDEX [17] experiments using Ge as targets, as as 
well as the KIMS experiment [T8] using CsI(Tl) crystals, 
are in disagreement with some or all of these claims. 

To achieve sensitivities to WIMPs beyond the current 
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experimental bounds, detectors with larger targets, lower 
background, and lower energy threshold are required. In 
the past decade, dual phase xenon detectors have rapidly 
emerged as one of the most promising technologies in 
WIMP direct detection, leading the WIMP search sen¬ 
sitivity in a wide range of parameter space mmm, 
demonstrating superior scalability in mass, and the capa¬ 
bility to shield against and reject background. However, 
in comparison to the cryogenic bolometers Bang or semi¬ 
conductor ionization detectors m eh nzi, dual phase 
liquid xenon detectors have not demonstrated the abil¬ 
ity to obtain a comparably low energy threshold. Con¬ 
ventionally, the issue is attributed to insufficient light 
collection efficiency or a lack of understanding of the 
low energy nuclear recoil (NR) quenching factor. In re¬ 
cent years, the LUX and PandaX collaborations oper¬ 
ated newly designed liquid xenon detectors which were 
constructed to optimize light collection efficiency. At the 
same time, a comprehensive model of scintillation and 
ionization processes in xenon known as the NEST EOT 
[22] , developed with simple phenomenological models and 
based on consideration of world data, is gradually be¬ 
ing adopted in the xenon field. The values of the rela¬ 
tive scintillation efficiency (L e ff) from NEST decrease 
continuously down to zero energy, which is consistent 
and slightly lower than that from an independent phe¬ 
nomenological calculation [23^. These developments call 
for careful re-examination of the low mass WIMP sensi¬ 
tivity using xenon detectors. In Ref. m, we reported 
the first 17.4 live days null search results in PandaX-I. 
In this paper, we present an improved analysis including 
the full PandaX-I data set, starting from May 26, 2014 
to Oct. 16, 2014, with a total of 80.1 live-day exposure 
in the search for dark matter. We shall refer to these 
data as dark matter search data in the remainder of this 
paper. 

THE PANDAX-I EXPERIMENT 

PandaX is a dual-phase liquid xenon dark matter ex¬ 
periment [[24] located at the China Jin-Ping Underground 
Laboratory (CJPL) [25]. The first phase PandaX-I is a 
pancake-shaped 120 kg detector optimized for light col¬ 
lection targeted for low mass WIMPs ns mi- The xenon 
chamber is a stainless steel inner vessel with an inner di¬ 
ameter of 750 mm, housing approximately 450 kg of liquid 
xenon. The entire inner vessel sits in an outer vacuum 
vessel constructed from 5-cm thick high-purity oxygen- 
free copper serving also as a radon barrier and electro¬ 
magnetic shield, and enclosed by a passive shield made of 
copper, polyethylene, lead, and polyethylene, from inner 
to outer layers. The gap between the outer vessel and the 
passive shield is continuously flushed with boil-off nitro¬ 
gen to maintain a radon level of less than 5 Bq/m 3 , more 
than a factor of 20 below the level in the experimental 


hall. The central time-projection-chamber (TPC) is a 
cylinder with diameter of 60 cm and height of 15 cm con¬ 
fined by a cathode grid (—15 kV) at the bottom, a gate 
grid (—5kV) and an anode mesh (ground) separated by 
8 mm, below and above the liquid level respectively, and 
a surrounding Polytetrafluoroethylene (PTFE) reflective 
wall. After a particle-xenon interaction, prompt scin¬ 
tillation photons (SI signal) are produced in the liquid. 
Ionized electrons are then drifted vertically upward by 
an induced drift field and extracted into the gas by an 
extraction field, producing the electroluminescence (S2 
signal). A top photomultiplier tube (PMT) array con¬ 
sists of 143 Hamamatsu R8520-406 (1-in square) tubes 
and a bottom array holds 37 Hamamatsu R11410-MOD 
(3-in circular) tubes. The PMTs view the active volume, 
collecting photons from both the SI and S2 signals, with 
the bottom array dominating the light collection for both 

51 and S2 signals. The radioactivity from the bottom 
PMT array is shielded by a layer of 5 cm thick LXe be¬ 
low the cathode and the PMT surface. The average dark 
rate per tube, i.e. the rate of random single photoelec¬ 
trons (PEs), is 0.06 kHz and 1.07kHz for top and bottom 
PMTs, respectively. The time separation between SI to 

52 signals gives the vertical position of the interaction, 
and the horizontal position is encoded in the S2 charge 
pattern in the PMT arrays. Multiple scatter events can 
be identified from the data by events which contain mul¬ 
tiple S2 signals, either separated in time if they happen 
at different vertical position or separated in the horizon¬ 
tal plane if there are multiple charge clusters in the PMT 
pattern. Gamma ray background produces electron re¬ 
coil (ER) events whereas the dark matter signal produces 
nuclear recoil events. The ratio of SI and S2 signal area 
gives a powerful means of ER rejection when looking for 
DM-like NR signals [26] . 

The PMT waveforms, amplified by a factor of 10 us¬ 
ing Phillips 779 amplifiers, are recorded by CAEN V1724 
14-bit 100 MS/s digitizers. The trigger for the data ac¬ 
quisition system (DAQ) is generated based on the major¬ 
ity outputs from the five digitizer boards for the bottom 
PMT arrays. For low energy signals in the dark mat¬ 
ter region, the trigger is generated by S2 with a thresh¬ 
old of about 89 total PE, whereas higher energy events 
were triggered primarily by SI with a charge threshold of 
about 65 PE. Each readout window is 200 /us long, with 
approximately equal division of pre- and post-trigger 
readout times. The PMTs are balanced to a gain of 
2 x 10 6 , with a recorded amplitude of the single photoelec¬ 
trons roughly at 60 digitizer bits. To save data volume, 
segments with waveform samples less than 20 digitizer 
bits from a pre-loaded baseline are zero-suppressed. For 
non-suppressed segments, 40 time samples before and af¬ 
ter the 20 bit threshold crossing are recorded. 

Three types of data runs were taken during the 
PandaX-I running period, the WIMP search, 60 Co ER 
calibration, and 252 Cf NR calibration runs. A summary 
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of the data taken is given in Table |T] Various cuts (dis¬ 
cussed below) are applied to remove periods with unsta¬ 
ble operating conditions, leading to a difference between 
DAQ time and the live time. 


Run type 

DAQ Time 
(hr) 

Live Time 
(hr) 

Trigger Rate 
(Hz) 

DM 

2,158.32 

1,923.11 

3.58 

252 Cf 

95.32 

94.05 

17.95 

60 Co 

405.14 

361.47 

22.23 


TABLE I: Summary of data taken during the entire PandaX-I 
running period. 

Two independent analyses were developed within the 
collaboration, utilizing different signal window selection, 
signal identification and reconstruction, event selection 
cuts and efficiencies, as well as the final fitting method. 
The two analyses were thoroughly cross checked at var¬ 
ious analysis stages, yielding consistent results. In the 
remainder of this paper, we will elaborate one of the anal¬ 
yses, and the other one is detailed in Ref. m- 

DATA PROCESSING AND SELECTION CUTS 

A number of improvements have been made in the data 
analysis pipeline compared with the first results [T3]. We 
shall describe the general procedure in steps below, with 
major improvements highlighted. 

The raw data files are screened for basic data quality 
before being processed for physics analysis. Detection of 
PMT high voltage outages is applied to filter data sets 
with low light collection. A nominal trigger rate below 
10 Hz is required to reject those data sets which are seri¬ 
ously contaminated by noise during times when running 
conditions are poor. Files with unexpected discharges 
from electrodes can be discriminated using the average 
number of S 1-like and S2-like signals in a waveform. If 
containing an average of larger than 40 Sl-like or 10 82- 
like signals, the events will be removed. Dark rates from 
PMTs are tracked and used to characterize the stabil¬ 
ity of the detector. A low random coincidence rate is 
essential and a cut is developed on PMT dark rates to 
minimize contamination. 

Baseline subtraction is performed on each waveform. 
In this analysis, the baseline is calculated based on the 
pre-samples from each waveform segment to suppress the 
drift and overshoot of baselines, whereas in Ref. m only 
weekly calibrated baselines were loaded. This update 
caused a downward shift of the light yield of approxi¬ 
mately 6% at 40keV ee electron-equivalent energy. 

Several malfunctioning PMTs are inhibited in the anal¬ 
ysis. During operation, four bottom PMT channels grad¬ 
ually developed connection problems, manifest as im¬ 
proper base resistance or capacitance, and were inhib¬ 
ited in the analysis to avoid a time dependent light yield. 


Among the rest of the bottom PMTs, a number of them 
experienced excessive dark rate (10 kHz and above) dur¬ 
ing the run but could sometimes be recovered through 
power cycling or lowering the corresponding high volt¬ 
age. One channel was fully inhibited due to unstable 
dark rate. The channel inhibition led to another 10% re¬ 
duction in light yield. On average, two to three bottom 
PMTs had to run at a lower gain (< 1 x 10 6 ) to main¬ 
tain a manageable dark rate. For the top array, seven 
PMTs gradually developed problems during the run and 
were inhibited as the problem showed up, but has lesser 
impact to the analysis presented here. 

Gain correction was applied to baseline-subtracted 
waveforms based on the the results of weekly LED cali¬ 
bration runs. A hit finder algorithm identifies signal hits 
channel-by-channel while tagging noise primarily due to 
the periodic 200 kHz electromagnetic interference from 
the CAEN PMT high voltage supplies, occasionally fluc¬ 
tuated above the 20 bit zero suppression threshold. The 
waveform is then integrated in the “hit window” to de¬ 
fine the hit charge. The hits in different channels are 
then clustered in time using an improved charge depen¬ 
dent algorithm with high efficiency for in-time short SI 
signals while avoiding splitting a low charge but wide S2 
into multiple clusters. For each cluster, a software sum 
is formed on all digitized channels, from which one com¬ 
putes the full-width-half-maximum, the full-width-1/ 10- 
maximum, as well as the number of peaks in the cluster. 
A binary decision tree method is developed to sort any 
given cluster into Sl-like and S2-like signals based on 
these variables. The identification efficiency is verified 
to be nearly 100% by checking waveforms of thousands 
identified clusters by eye. 

We developed further signal-level cuts to identify spu¬ 
rious noise in the Sl-like and S2-like signals. For Sl-like 
signals, we employed a further ripple-pattern cut on the 
software summed waveform, a cut on the ratio of charge 
computed from the summed waveform to the total hit 
charge , and a cut on the ratio of the total number of 
noise hits to the total hits in the cluster. In addition, 
cuts are placed on the ratio of the height to area and 
that of the height to width. To avoid signals from af¬ 
terpulsing, SI signals are required to be before the first 
good S2 signal. For S2, we developed a shape symmetry 
cut to remove events very close to the anode when SI and 
S2 cannot be easily separated in time, identifying those 
events with a characteristic sharp spike at the beginning 
of the S2 pulse. In addition, S2-like signals will be dis¬ 
carded if its ratio to the largest S2 is less than 1% or if 
such signal is consistent with a single-electron S2 with 


1 The summed waveforms for the 200 kHz noise tend to show a 
clear ripple feature, leading to a cancellation in the corresponding 
charge. 
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charge less than 30 PE; the inefficiency due to these two 
cuts is estimated to be negligible by a Monte Carlo (MC) 
simulation. 

Periods with unstable overall PMT signal rates were 
identified during the run, possibly due to small discharges 
on the electrodes, producing light pulses. In addition to 
the earlier file-by-file cuts, a number of tighter data qual¬ 
ity cuts were placed event-by-event to search for “dirty 
waveforms”, including a cut on the total number of Si- 
like signals, the ratio of SI to total charge before the 
first S2, and the ratio of the sum of SI and S2 to the 
total charge (which also suppresses multiple scattering 
events). To avoid ambiguity due to multiple S 1-like sig¬ 
nals per event, we require that the number of SI should 
be either 1 or 2, and in case of the latter the maximum 
SI is identified as the primary if the charge of the second 
SI does not exceed 50% of the first. Finally, to suppress 
accidental background, we placed a 3-hit coincidence cut 
on any good signal, and a 300 PE cut on S2 (discussed 
later). 

Finally, SI and S2 signals are reconstructed into physi¬ 
cal events. The vertical position is determined by the sep¬ 
aration between SI and S2 signals, assuming a drift speed 
of 1.7mm//xs under the drift field of 0.67kV/cm IT51I28] . 
Horizontal position of the interaction is reconstructed 
with the S2 PMT charge pattern using multiple algo¬ 
rithms. As an improvement to the charge center of grav¬ 
ity (CoG) method, we developed a fast charge pattern 
template matching (TM) method. The expected charge 
templates were generated using a custom GEANT4 [29] 
based Monte Carlo which simulates optical photon prop¬ 
agation from S2 signals in the PandaX-I TPC geometry, 
identical to the templates used in the fast artificial neu¬ 
ral network (FANN) reconstruction method developed by 
the independent analysis. The difference in the horizon¬ 
tal positions from the FANN and TM methods is on av¬ 
erage 5 mm, obtained using 40 keV de-excitation events 
from neutron calibration data. This is independent of the 
radial and vertical positions and consistent with the ex¬ 
pected position resolution from MC. In the analysis pre¬ 
sented here, the TM reconstruction is chosen. To identify 
multiple scattering at the same vertical location, we set a 
charge clustering cut by requiring the horizontal distance 
between the CoG and TM positions be less than 55 mm 
apart. 

DETECTOR CALIBRATIONS 

The PandaX-I detector has been carefully calibrated 
using various methods to perform an effective search for 
low-mass WIMPs. Single electron events were identified 
to calibrate the single electron gain (SEG) in the elec¬ 
troluminescence. Neutron-activated X-rays were used to 
determine the photon detection efficiency for SI (PDE) 
of the PMTs which signifies the sensitivity of our detec¬ 


tor in the low-mass region, and the electron extraction 
efficiency (EEE) from the liquid. A neutron calibration 
with 252 Cf was used to generate the NR events that were 
used to define the DM search window. Finally, a gamma 
calibration with 60 Co was used to find the leakage of ER 
background into the search window. 

Neutron calibrations have been taken several times 
throughout the run with a total exposure of 95 hrs (see 
Table I). The 40keV ( 129 Xe) and 80keV ( 131 Xe) inelas¬ 
tic recoil X-rays are used in calibrating the uniformity 
for both SI and S2 in the detector. For a fixed energy 
deposition in the detector, the PMT arrays see different 
light and charge yields depending on the spatial loca¬ 
tion of the event, which must be corrected to a detector 
average before further analysis is performed. The uni¬ 
formities for SI and S2 are verified to be decoupled in 
the vertical direction and horizontal plane. The horizon¬ 
tal variation of 40keV S2 peaks in the 54.0 kg fiducial 
volume is measured to be ±36%, which dominates the 
detector non-uniformity. The vertical uniformity for the 
S2 signals, characterized by an exponential “electron life¬ 
time”, reflects the electronegative impurity level in the 
detector which tends to attenuate the charge signal dur¬ 
ing drifting. The average electron lifetime in our detec¬ 
tor is fit with a decaying exponential and determined to 
be 328±8 /is (an attenuation length of about 60 cm as 
compared to the 15cm maximum drift distance). On 
the other hand, the variation of the SI peak in the fidu¬ 
cial volume in the vertical (horizontal) direction is ±8.5% 
(±9.5%). All discussions in the remainder of the paper 
are made with the uniformity corrections taken into ac¬ 
count. 

One of the important properties of the detector is 
SEG, the average number of PEs observed in PMTs from 
single-electron electroluminescence. It can be determined 
from the PE distribution of the smallest S2 signals taken 
at any normal detector run, fit with a double Gaussian 
function with means related by a factor of 2 from the 
charge quantization (see Fig. [I]). The SEG is determined 
to be 18.4±1.6 PE/e, where the uncertainty is estimated 
by varying parameters in the hit clustering algorithm as 
well as the fitting function and range. 

The PDE and EEE can be determined from the 40 
and 80keV X-ray events during the neutron calibration. 
The events collected are shown in the S2 vs. SI plot in 
Fig. i The location of the 40 keV peak (with decay time 
less than Ins) is at 178.8 PE in SI, with an average 11.6 
PE mixture from the associated NR, estimated from the 
pure NR events seen at low energy as well as through 
MC. At this energy, our detector has a SI photon yield 
of 4.2 PE/keV. Using the NEST-0.98 model [21], this 
corresponds to 6.0 PE/keV at zero electric field at the 
standard 122 keV, in comparison to the 3.9 PE/keV ob¬ 
tained XENON100 [IT] and 8.8 PE/keV in LUX p2]. 

The electron-equivalent energy of a given events can 
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FIG. 1: The PE distribution of the smallest S2 signals (cor¬ 
rected for horizontal non-uniformity), summed over all top 
and bottom PMTs. The single electron gain is determined by 
fitting two constrained Gaussians, shown as the dashed green 
lines. 



SI [PE] 


FIG. 2: Distribution in uniformity-corrected SI and S2 for the 
de-excitation gamma peaks in the neutron calibration runs. 
The anti-diagonal lines are the anti-correlation fit at the two 
energies. The dashed vertical and horizontal lines indicate 
the mean NR energy in mixture with the gamma energies at 
40 and 80 keV. 

be reconstructed from the light and charge outputs as 

E rec = (Sl/PDE + S2/SEG/EEE) x W, (1) 

where E rec is the reconstructed energy in keV ee splitting 
into scintillation and ionization parts, and W = 13.7 eV 
is the average energy to produce a scintillation photon or 
to liberate an electron m- The anti-correlated fluctua¬ 
tions in the light and charge outputs due to electron-ion 
recombination is naturally accounted for in Eq. [l] Sim¬ 


ilar to Ref. na. we performed anti-correlation fits using 
Eq.0 to the 40 and 80keV de-excitation peaks, as well 
as the neutron-induced meta-stable 129m Xe (164 keV) de¬ 
cay gamma rays after the neutron calibrations 2 . The 
PDE (EEE) determined with the 40keV ee peak is 9.6% 
(82.1%). The fractional uncertainties are estimated to 
be 10% and 9%, respectively, based on the difference in 
values obtained at other two energies, as well as those in 
Ref. [151 . 

To facilitate the comparison of our data with model 
prediction, we convert the peaks in SI and S2 into a per 
unit energy total photon yield (L y ) and charge yield (C y ), 
using 

L y = (51)/PDE/£ rec , 

C y = (S'2)/SEG/EEE/E rec , (2) 

where (51) and (52) here refer to the location of corre¬ 
sponding peaks in the distribution. In Fig. j3j our mea¬ 
sured data is compared to the mean values in NEST- 
0.98 [5Tj under the same drift field. Reasonable agree¬ 
ment is found at all four energy peaks in 252 Cf data (40, 
80, 164, 236 keV). The uncertainties shown in the fig¬ 
ure, aside from the statistical uncertainties in the peak 
determinations, arise from the systematic uncertainties 
of the PDE and EEE determination through the anti¬ 
correlation fits. 



FIG. 3: The L y (blue) and C v (red) (in units of quanta per 
keVee) extracted based on PDE and EEE obtained at dif¬ 
ferent energies overlaid with corresponding curves predicted 
by NEST-0.98. The reconstructed energy spectra for the de¬ 
excitation peaks and meta-stable xenon isotopes are also over¬ 
laid with y axis scaled for visual clarity with fitted energy 
resolutions indicated in the figure. 

In the 252 Cf NR calibration runs, the single events at 
very low energy with Sl< 30 PE are expected to have 


2 We did not perform anti-correlation fits for the 131m Xe 236 keV 
gamma lines since it was difficult to separate the peak cleanly 
from the background. 
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less than 1% contamination from the ER band based 
on MC simulations, and the latter can therefore be ne¬ 
glected. The distribution of these low energy events in 
logio(S2/Sl) vs. SI is shown in Fig. |4ja). It can be seen 


NEST-1.0 model [22] is used instead, the MC can also be 
tuned to agree with the data by increasing the EEE up by 
2%, much less than its assigned uncertainty. The tuned 
MC is used as the true physical distribution to extract 
the NR efficiency. 



SI x,y,z corrected [PE] 


FIG. 4: The band of logio(S2/Sl) versus SI for the NR cali¬ 
bration data without (a) and with (b) the “X” cut. See text 
for definition of the cut. The 54.0-kg fiducial cut is applied. 
The solid blue line in (b) is the median of the pure NR band 
in MC, and the dots are the Gaussian mean obtained from 
the data for Sl> 15 (where the detection efficiency is flat so 
that the data and MC comparison is straightforward). The 
dashed magenta lines in both figures are the 300 PE cut on 
S2, below which no dark matter candidate is considered. The 
gray dashed lines are the equal energy lines with NR energy 
indicated in the figures. 

that there are scattered events with suppressed S2, pro¬ 
ducing an asymmetric NR band. Based on the charge 
pattern of SI signals, it was determined that such events 
(called “X” events [3D]) are due to multiple scattering 
of neutrons with some energy deposition in the “charge¬ 
less” region, either below the cathode, or in the xenon 
“skin” between the PTFE wall and the stainless steel in¬ 
ner vessel (viewed partially by the outermost ring of the 
top PMT array). They have to be properly taken into 
account to correctly calibrate the NR efficiency. 

To compare the data with expectation, a GEANT4 MC 
is developed to simulate the 252 Cf runs, which produces 
both single-scatter pure NR and “neutron-X” recoil spec¬ 
tra, and employs the NEST-0.98 nuclear recoil model m 
with the PDE, EEE and the SEG obtained above. After 
global tuning of the strength of the ” neutron-X” events 
in the MC, excellent agreement is found in different slices 
of SI between the data and MC (Fig. [ 5 ]) 3 . If the new 


2 to 5 PE 5 to 10 PE 



10 to 15 PE 



FIG. 5: Comparison of the distribution of log 10 (S2/Sl) in 
the NR data and tuned MC in six slices of SI as indicated 
by the figure titles: data (blue), tuned “neutron-X” events 
in the MC (magenta), and the sum of pure NR and tuned 
“neutron-X” in MC (red). In each slice of SI, the value of 
log 10 (S2/Sl) is shifted relative to the median value in that 
slice. The efficiency in Fig. [6] has been applied to the MC to 
compare with the data. 


To suppress the “X” events, a charge asymmetry cut 
between the top and bottom PMT arrays as well as a 
cut on the ratio of the maximum single PMT charge to 
the total on SI were applied to all data including 252 Cf, 
60 Co, and DM data sets. The NR distribution after the 
cut is shown in Fig. [ 4 ) 3 , where the low S2 “X” events are 
significantly reduced. Our overall analysis cut efficiency 
for NR events with SI > 10 PE is estimated by comparing 
the number of 252 Cf events in (S1,S2) bins before and 
after all cuts in this energy region, and an approximately 
uniform 77.5% value is obtained. At lower energy, the 
overall NR efficiency is estimated by taking the ratio of 
the measured distribution to the tuned MC, anchored at 
77.5% at higher energy. The resulting two-dimensional 
distribution of the NR efficiency is shown in Fig. [6] 

The ER calibration is performed with a 60 Co gamma 
source, interleaved frequently during dark matter data 
taking. Low-energy 7 rays are produced through the 
well-known Compton scattering mechanism. The distri¬ 
bution of the single scatter ER events in logio(S2/Sl) vs 
SI is shown in Fig. 7(a) All cuts, including a fiducial 


3 A fluctuation of 17% from the gas gain, in addition to the nomi¬ 


nal statistical fluctuations introduced by NEST, helps to match 
the measured width in each SI slice. 
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FIG. 6: The nuclear recoil efficiency in (S1,S2) determined 
using the method discussed in the text. 


cut, have been applied. For events with SI between 2 
and 30 PE, 12 out of 1,520 events were located below 
the median of NR in logi 0 (S2/Sl). Subtracting the ex¬ 
pected 1.65 events from accidental coincidence (see later 
discussions), the remaining ER leakage is 0.68±0.23% of 
the total, consistent with a pure Gaussian expectation 
(0.5%) obtained by fitting the ER band distribution. 

Within the SI range of 10 to 30 PE, the efficiency for 
ER detection and selection is estimated to be 75.7% by 
taking the ratio between the final number of events after 
all cuts and the raw events on the ER band. At lower 
energy, the efficiency is estimated by taking the ratio 
between the measured and expected SI spectrum from 
MC with 75.7% at higher energy as an anchor (shown in 
Fig. |7(b)| . The overall efficiency is approximately 71.5% 
in entire 2-30 PE range. 


BACKGROUNDS IN DARK MATTER SEARCH 
DATA 


The low-energy dark matter window was blinded in the 
analysis until all data cuts were determined. The cuts on 
SI and on the fiducial volume were optimized from a 
figure-of-merit based on the expected below-NR-median 
backgrounds of the ER, the accidental background (sta¬ 
tistically determined from data), and the neutron back¬ 
ground (MC estimates). The final optimized search win¬ 
dow on SI is from 2 to 30 PE, and that on S2 is 300 to 
10,000 PE. The fiducial cut is determined as r 2 < 500 cm 2 
with a drift time between 10 to 80 fx s, resulting in a fidu¬ 
cial mass of 54.0±2.3kg. In what follows, we shall discuss 
the background contributions in the dark matter search. 

ER background Expected ER background in our fi¬ 
nal candidate sample with all cuts imposed, summarized 
in Table [TIJ has been estimated with a GEANT4-based 




(b)ER efficiency 


FIG. 7: (a) The band of log 10 (S2/Sl) versus SI for the 60 Co 
ER calibration runs with the median and ±2cr of the band 
indicated as the solid and dashed blue lines, respectively. The 
median of the NR band is indicated as the solid red line, below 
which the 12 leaked events are plotted as green markers. The 
dashed magenta line is the 300 PE cut on S2. The gray dashed 
lines are the equal energy lines with ER energy indicated in 
the figures; (b) The ER efficiency in SI obtained by taking 
the ratio between the data and MC (histogram), and the red 
curve is a fit to the efficiency. 


MC program, with a few updates compared to that in 
Ref. JT3]. First, by taking into account the additional en¬ 
ergy deposition in the below-cathode region (“X” events) 
observable through PMT arrays, some of the MC events 
shifted above of the dark matter search window, leading 
to a reduction of background from almost all components. 
Second, the radioactivity level of the stainless steel vessel 
was updated with a counting measurement with much 
better statistics, also resulting in a reduction in back¬ 
ground expectation. Third, the internal 85 Kr, 222 Rn and 









































Source 

background level (mDRU) 

Top PMT array 

4.7±2.3 

Bottom PMT array 

2.3 ±1.5 

Inner vessel components 

3.8 ±2.2 

TPC components 

1.9 ±0.9 

85 Kr 

2.6 ± 1.2 

222 Rn & 220 Rn 

0.5 ± 0.2 

Outer vessel 

0.9±0.6 

Total expected 

16.7±3.9 

Total observed 

23.6±3.5 


TABLE II: The expected and observed background rates in 
the fiducial volume and in the dark matter search window. 
mDRU = 10 -8 evt/day/kg/keV ee . Uncertainties in the MC 
prediction originate from uncertainties in the material ra¬ 
dioactivity screening, except those for Rn and Kr which are 
due to the uncertainties in the PandaX data. 


220 Rn levels were studied with the statistics of the full 
dark matter search data sample with the same delayed 
coincidence techniques as in Ref. m The measured Kr 
concentration in Xe is 68 ±29 ppt mole/mole (uncertain¬ 
ties mainly due to event selection methods in the anal¬ 
ysis) assuming a 2 x KR 11 isotopic abundance of 85 Kr, 
leading to an expected background of 2.6 ± 1.2 rnDRU 
based on the MC. The 222 Rn and 220 Rn backgrounds 
were determined to be 0.7 ± 0.2 and 0.15 ± 0.06 mBq 
in the fiducial volume, respectively, with uncertainties 
primarily arising from event selection cuts. The result¬ 
ing background of 0.5 ± 0.2 mDRU in the dark matter 
energy and fiducial volume search window is estimated 
by MC with an improved treatment taking into account 
the non-secular equilibrium due to the long-lived isotope 
210 Pb (r = 22.2 year). The overall ER background in 
the dark matter search data estimated from radioactiv¬ 
ity counting is 16.7±3.9 mDRU. This is consistent with 
the ER background 23.6 mDRU extrapolated from events 
with Sl> 30 after efficiency correction (±15% depending 
on the energy cut as well as efficiency modeling), assum¬ 
ing a flat distribution of the ER background in keV ee at 
very low ER energy based on the MC. 

The real relevant ER background for dark matter 
searches is formed from events that leak below the NR- 
median, including those due to detector effects as well 
as the so-called “gamma-X” events with partial energy 
deposition in the “chargeless” regions. To reliably esti¬ 
mate the number of such events, it is best to use the ER 
calibration data where such events are included with the 
right proportion. 

Neutron background The neutron background is esti¬ 
mated using a combination of SOURCES-4A [3T, and 
GEANT4 simulation, leading to an estimate of 1.45 
events within the 54.0x80.1 kg-day exposure before effi¬ 
ciency cuts, and about 0.35 events after all cuts. This 
yields 0.18 neutron background events below the NR 
medium line. We assign a generous 50% uncertainty to 
the MC estimate. Alternatively, a 90%-confidence-level 


upper limit of 1.15 neutron events can be set based on 
the single to multiple NR scattering ratio from the MC 
and the absence of the multiple scattering NR’s in the 
dark matter search. 

Accidental background In our dark matter search 
data, we find a significant number of isolated SI and S2 
events, which yield a substantial background. An isolated 

51 is an event occurring without an obvious S2 nearby. 
These signals are likely from multiple origins, e.g. light 
leaking into the TPC due to interaction in the skin re¬ 
gion, small discharges in the TPC due to impurities or 
high voltage, and the accidental coincidence of SPE be¬ 
tween PMTs. An isolated S2 is an event without an SI 
proceeding the waveform, which can be due to events 
with very low energy of which SI cannot be detected. 
In addition, based on a visual inspection of isolated S2 
events, it was noticed that a significant fraction of such 
events have a spiky timing profile at the beginning of 

52 (but cannot be efficiently rejected with existing al¬ 
gorithms), implying that these S2 events happened very 
close to the gate grid where SI and S2 can no longer be 
separated. 

In our dark matter data, isolated SI events are es¬ 
timated by looking for uncorrelated SI events before a 
large SI (which is associated with a trigger) yielding a 
rate of about 23 Hz, with the charge distribution shown 
in Fig. |8(a)| Isolated S2 events are measured with a rate 
of about 240 events/day for S2 within 300 to 10,000 PE 
(the 300 PE cut is imposed to balance between the sup¬ 
pression of such background and the loss of low-energy 
sensitivity), without obvious non-uniformity in the hori¬ 
zontal plane. The charge spectrum of such S2 signals is 
shown in Fig. 8(b) The rate of such events in the 60 Co 


calibration runs increased to about 568 events/day, and 
the amount of the rate increase is consistent with the MC 
expectation. 

Isolated SI and S2 events produce accidental coinci¬ 
dences which mimic real events. Such background can 
be statistically evaluated by forming random pairs in SI 
and S2, and the resulting distribution in logio(S2/Sl) 


versus SI is shown in Fig. 8(c) The overall rate in the 


dark matter data is estimated to be 35.1 events in 80.1 
day with a conservative 10% uncertainty based on the 
statistical uncertainty of a day-long dark matter search 


CANDIDATE EVENTS FROM 80.1 DAY DARK 
MATTER SEARCH DATA 

For the dark matter search data, the event rates after 
different levels of cuts are summarized in Table HU The 
data quality cuts remove a large fraction of the multiple 
scattering events, reducing the total number consider¬ 
ably, which also explains that the subsequent single-site 
cut has a small effect on the remaining number of events. 
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SI [PE] 

(c) Accidental 


FIG. 8: The measured isolated SI (a) and S2 (b) charge 
spectra, and the band of log 10 (S2/Sl) versus SI (c) for the 
statistically obtained accidental background, with the sharp 
cutoffs at the top and bottom corresponding to the 10,000 
and 300 PE cut. The median of the NR band is indicated as 
the red line. 
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Within prescribed cuts, 542 events were found in 54.0 kg 
x 80.1 days. The event distribution in r 2 vs. drift time in 
the TPC is shown in Fig. 9(a) The event projections in 
r 2 (with two position reconstruction methods) and drift 
time are also compared to the expected ER distribution 
from the Monte Carlo, where good agreement is achieved. 


FIG. 9: (a) The vertex distribution of events in the TPC dur¬ 
ing dark matter data taking. The blue dashed box indicates 
the fiducial volume, and the red dashed box indicates the en¬ 
tire active volume within the TPC, confined by the cathode, 
gate grid, and the PTFE wall. The events below the NR me¬ 
dian are indicated by the green markers; (b) projection in r 2 , 
data (TM method for position reconstruction) in blue, data 
(CoG method) in green, and the MC energy deposition CoG 
in red; (c) projection in the drift time, data (TM) in blue and 
the MC energy deposition CoG in red. 
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Cut 

# Events 

Rate (Hz) 

All triggers 

24,762,972 

3.58 

Quality cut 

6,127,280 

0.88 

Single-site cut 

5,050,845 

0.73 

SI range (2-30 PE) 

62,872 

9.08 x 10~ 3 

S2 range (300-10,000 PE) 

44,171 

6.38 x 10“ 3 

Fiducial volume 

542 

7.83 x 10~ 5 


TABLE III: The event rate in the DM run after various level 
of analysis cuts. 


The distribution of events in logio(S2/Sl) vs. SI is 
shown in Fig. 10 The majority of the events are consis¬ 
tent with an ER origin. The events located higher than 
the ER band at low SI are the accidental backgrounds, 
more prominent than those in the 60 Co calibration run 
due to the much lower ER event rate in the dark mat¬ 
ter search data. Seven of the candidate events are lo¬ 
cated below the median of the NR band indicated by 
the green markers in Figs. 9(a) and 10 For comparison, 
the expected background in the total sample as well as 
those below the NR mean is shown in Table IIVI The ER 
background is estimated based on the 23.6 mDRU value, 
a corresponding ER energy range of 6.9keV ee , and an 
average ER efficiency of 71.5%. The below-NR-median 
accidental background is estimated based on the distri¬ 
bution in Fig. 8(c) Summing over all the contributions, 
we expect 6.9 background events below the NR median. 
No significant excess above the background is observed. 



ER 

Accidental 

Neutron 

Total 

expected 

Total 

observed 

All 

503.7 

35.1 

0.35 

539.1 

542 

Below 
NR med 

2.5 

4.2 

0.18 

6.9 

7 


TABLE IV: The expected and measured events (in units of 
events) in 80.1 live-day dark matter search data. 


FITTING METHOD 


To maximally use the information from the data, in¬ 
stead of choosing only the below-NR-median region to 
search for DM like in Ref. [T3j , in this analysis we defined 
a much extended DM window with SI between 2 and 30 
PE and S2 between 300 to 10,000 PE. To fit all data, an 



FIG. 10: The band of log 10 (S2/Sl) versus SI for the dark 
matter search data. The ER band is indicated by the blue 
solid line (median) and the dashed line (±2cr). The median 
of the NR band is indicated as the solid red line. The dashed 
magenta line is the 300 PE cut on S2. The green stars rep¬ 
resent events below the NR median. The gray dashed lines 
are the equal energy lines with NR energy indicated in the 
figures. 


unbinned extended likelihood function is constructed as 
C = Poisson(N m |N exp ) x 

TTi = Nm <N dm (1 + S dm )Pdm( Sl\ S2‘)e^(Sl\ S2*) 

111=1 ^ N 

1 y exp 

| AWl + fefl)P.E fl (SP,S2*) 

N e xp 

| N Acc (l + 5 Acc )P Acc (Sl i ,S2*) 

N GX p 

N n bkg{ 1 + S n bkg)Pnbkg(Sl\ S2*)ejVij(Sl*, S2*) 

Nexp J 

x G{S dm , 0.2 )G(S E r, 0.15 )G(6 Acc , 0.1 )G(S nbkg , 0.5), 

(3) 

where N m and N exp are the total number of measured 
and fitted candidates with 


Nexp = N DM (eNR)DM( 1 + 5 dm) + N E r{ 1 + S E r ) 

T -N.4cc(1 + Sacc ) T N nb kg(cNR)nbkg{ 1 “b 5nbkg) • 

(4) 

As indicated in Fig. [9j the position dependence of events 
in the fiducial volume is rather weak and is therefore 
ignored here for simplicity. N E m ( N nb k g ) is the total 
number of WIMP particles (neutrons) interacting with 
the detector during the measurement before efficiency 
and acceptance cuts. N E m is computed for each given 
pair of WIMP mass and cross section (m x ,a n - x ) assum¬ 
ing the isothermal DM halo model [32l 331 with a local 
dark matter density of 0.3 GeV/c 2 /cm 3 , a circular veloc¬ 
ity of 220km/s, a galactic escape velocity of 544 km/s, 
and an average earth velocity of 245 km/s. Pdm( SI 4 , S2 1 ) 
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and P n bkg( SI 1 , S2 1 ) are the probability distribution func¬ 
tions (PDFs) of NR recoil signals for a WIMP with given 
mass and neutron background, respectively, obtained us¬ 
ing the NEST-based MC simulation employing the PDE, 
EEE, and SEG described earlier. ejvfl(Sl 4 , S2 l ) is the 
NR detection efficiency from Fig. [6j with acceptance set 
to zero if SI and S2 are outside the ranges of (2,30) 
and (300,10,000) PEs. To obtain the expected measured 
dark matter and neutron background events, the NR ef¬ 
ficiency function has to be averaged over the expected 
dark matter or neutron background PDF (( £ nr)dm and 
( £NR)nbkg ) in Eq. [ 4 ]). N er and N Acc are the total num¬ 
ber of ER and accidental background with detection ef¬ 
ficiency taken into account, and Per(S1 z ,S2 z ) (taken to 
be the same as that obtained from ER calibration from 
supported by the MC) and P J 4 CC (S1 I , S2 l ) (from 


7(a) 


8(c)) are the corresponding PDFs. The contamina- 


Fig. 

Fig. 

tion of the accidental background in the ER calibration 
run is neglected due to the dominating ER rate in the 
calibration runs. The expected background events are 
taken from the top row of Table |IV| To allow system¬ 
atic variation in the global efficiency, four normalization 
nuisance parameters (5dm, 5er, 5 Ac c and 5 n bk g ) are in¬ 
cluded for the four type of events, constrained by Gaus¬ 
sian variations (G’s in Eq. [ 3 ]) of 20% (DM), 15% (ER), 
10% (accidental) and 50% (neutron background) in the 
penalty terms |3J[ S3 • 

The average WIMP detection efficiency ( €nr ) dm , ob¬ 
tained by combining the NR efficiency with the WIMP 
PDF (with fluctuations in SI and S2 properly taken into 
account), strongly depends on the WIMP mass, which 
is depicted in Fig. El The lower the WIMP mass, 



FIG. 11: The overall detection efficiency for WIMP at differ¬ 
ent mass (red) and those obtained with PDE and EEE set at 
+lcr (solid green) and — la (dashed green). 

the softer the recoil energy distribution, therefore the 
selection threshold on SI and S2 would more strongly 
suppress the overall efficiency. To compare the effects 
of selection thresholds of different experiments on the 


NR energy, the mean NR energy curves in SI and S2 
are plotted in Fig. [12] based on the NEST-0.98 model 
with the PDE, EEE, and SEG values from PandaX, 
XENON100 [38], and LUX JT2]. Our selection thresh- 



FIG. 12: The comparison of the mean energy thresholds, 
translated from the cuts on SI and S2, for different exper¬ 
iments under the same energy model (NEST v0.98). The 
solid curves (red: PandaX, violet: XENON100, blue: LUX) 
represent the mean values for SI and S2 obtained from NEST 
with slanted ticks (along the equal-energy vector) indicating 
the corresponding mean NR energy in divisions of keV nr . The 
curve for PandaX based on NEXT vl.O is drawn as the green 
dashed curve. The selection thresholds in SI and S2 are in¬ 
dicated in the figure as the solid circles. The anti-diagonal 
dashed lines (red: PandaX, violet: XENON100, blue: LUX) 
are the equal-energy lines projected from the corresponding 
threshold points for different experiments. 

old is at about 4.2 k eV nr in both SI and S2. XENON100 
achieved a S2 threshold of less than 2 keV nr but a much 
higher SI threshold of about 8 keV nr . LUX, on the other 
hand, achieves an average 3keV rar . threshold on both SI 
and S2, but in Ref. [T2] they choose to drop the NR ef¬ 
ficiency entirely below 3keV nr . The NEST-1.0 model 
predicts a higher charge yield for NR, in which case our 
SI threshold would stay, but the S2 threshold would im¬ 
prove to about 2.8keV nr , leading to a better sensitivity 
for low mass WIMPs. Nevertheless, we chose the NEST- 
0.98 model to report our final WIMP results. 

The best fit value to maximize the likelihood func¬ 
tion is found at m x = 27.5 GeV/c 2 with a <T x ,jv at 
4.1xl0 -45 cm 2 . The value of the likelihood is also con¬ 
sistent with that from the null hypothesis within la, in¬ 
dicating no significant excess over the background. To 
set the WIMP search upper limit, a standard profile like¬ 
lihood ratio statistic is formed [33] [35]. A Feldman and 
Cousin approach [37] is used to fit the data as well as a 
large number of MC simulations using the signal hypoth¬ 
esis at each grid point of (m x ,cr n _ x ). The 90% c.l. up¬ 
per limit obtained with this approach is shown in Fig. |13| 
together with the world data, and is verified to be very 
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similar to that obtained assuming an approximate half-x 2 
distribution of the test statistic [34]. A binned likelihood 
method developed in the independent analysis yields an 
upper limit in good agreement with the above. The up- 



FIG. 13: The 90% c.l. upper limit for spin-independent 
isoscalar WIMP-nucleon cross section for the PandaX-I ex¬ 
periment (red curves). Recent world results are plotted for 
comparison: XENONIOO 225 day results QT] (black solid), 
LUX first results [12] (blue), SuperCDMS results [15] (or¬ 
ange solid), DarkSide results [T4J (magenta solid), CRESST- 
II 2014 limits [8] (brown dashed), and CDEX 2014 limits [lj] 
(solid violet). The claimed WIMP signals are shown as closed 
contours: CoGeNT 2014 results [6] (cyan solid), CDMS-II- 
Si results [5j (gold dashed), DAMA/LIBRA 3er contours |38] 
(green solid), and CRESST-II 2012 results [ 7 ] (brown solid). 


per limit excludes a WIMP mass of 10GeV/c 2 down to 
a cross section of 1.41 x 1CU 43 cm 2 , and the lowest ex¬ 
cluded cross section is 1.01xl0 _44 cm 2 at a WIMP mass 
of 44.7GeV/c 2 . Under the elastic, spin-independent, and 
isospin conserving WIMP-nucleon scattering model, our 
limits strongly disfavor the WIMP interpretation of the 
results from DAMA/LIBRA, CoGeNT, CDMS-II-Si and 
CRESST-II. It is noteworthy that the PDE and EEE 
used in this analysis are conservative in nature since 
we inhibited the unstable PMTs. In addition, we have 
considered the average WIMP detection efficiency with 
WIMP mass dependence in this analysis (Fig. 11). Com¬ 
pared to that in Ref. m , in which the DM efficiency is 
treated only as a function of SI, this treatment is more 
realistic. Even with these realistic treatments, our results 
still set a stringent limit at the low WIMP mass region, 
with a tighter bound than SuperCDMS above the WIMP 
mass of 7 GeV/c 2 , and the best reported bound in a dual 
phase xenon detector below a WIMP mass of 5.5 GeV/c 2 . 
Note that one of the key difference between this analysis 
and that from LUX in Ref. [12] is that the latter made 
a conservative choice to model no signal generation for 
events below 3 keV nr , while in our treatment the signals 
generation is continuous to zero energy therefore low en¬ 
ergy events below the mean energy threshold of 4.2 keV„ r 


could still fluctuate upwards into the detection region. 

The experimental sensitivity band is obtained using 
the same approach as above but with hundreds of 80.1- 
day background-only toy MCs based on Table [TV] using 
prescribed PDF for each event type, from which one ob¬ 
tains a distribution of “upper limits”. In Fig. 14 our 


upper limit is overlaid with the ±l-cr sensitivity band. 
Consistency is observed, confirming no significant excess 
over background. 

To study shape related systematic uncertainties sepa¬ 
rately 4 , we performed calculations of upper limits either 
by setting PDE and EEE both at +ler or — ler. The re¬ 
sulting limits are overlaid in Fig. [14] As expected, the 
higher efficiency would lead to tighter bounds in the low 
mass region and vice versa. The (more aggressive) upper 
limit obtained with dark matter PDFs generated from 
the NEST-1.0 model is very close to that with the +1 <t 
PDE/EEE. These are sizable influences but are compa¬ 
rable with the sensitivity band, therefore do not change 
the main conclusion of our results. 



FIG. 14: PandaX-I WIMP search limit from the data (red 
line) overlaid with the ±1 <t sensitivity band obtained from 
toy MC (yellow) as well as the alternative upper limits using 
either +lcr or —la values for the PDE and EEE, but with 
the same NEST-0.98 model. For comparison, a few world 
leading limits for the low mass WIMP are plotted: LUX 
first results [12] (blue), SuperCDMS results [15] (orange), and 
CRESST-II 2014 limits 0 (brown dashed). 


CONCLUSION AND OUTLOOK 


In summary, we report the low-energy dark matter 
search results with the 54.0x80.1 kg-day full exposure of 


4 The shape systematics could also be introduced into the fitter 
via nuisance parameters. However, to explicitly show the size of 
the effects and to simplify the fitter computation, we chose to 
apply these systematic variations “by hand”. 
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the PandaX-I experiment. In this analysis, compared to 
the first results, we made a number of improvements in 
signal identification, background classification and rate 
and shape estimates, a realistic treatment on the effi¬ 
ciency for very low recoil energy events, as well as profile 
likelihood ratio fits to obtain the final WIMP search limit. 
Observing no significant excess over background, our re¬ 
sults strongly disfavor the WIMP interpretation of the 
results from DAMA/LIBRA, CoGeNT, CDMS-II-Si and 
CRESST-II. Our bound is tighter than that from Super- 
CDMS above the WIMP mass of 7GeV/c 2 , and is the 
lowest reported limit below a WIMP mass of 5.5GeV/c 2 
in xenon dark matter experiments to date, showing that 
liquid xenon detectors can be competitive for low-mass 
WIMP searches. 

The results from PandaX-I are crucial in guiding the 
future development of the PandaX program. The sec¬ 
ond phase experiment, PandaX-II, constructed with a 
liquid xenon target of 500 kg sensitive mass and lower 
background materials for the cryostat and TPC, is un¬ 
der preparation at CJPL. The PandaX-II detector is ex¬ 
pected to improve both on the light and charge collection 
efficiency and push the dark matter sensitivity beyond 
the current best reach in a wide range of WIMP masses. 
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